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A wide spectrum of electrokinetic studies is modelled as isothermal ones to expedite analysis 
even when such conditions may be extremely difficult to realize in practice. As a clear and novel 
departure from this trend, we address the case of flow-induced electrohydrodynamics, commonly 
referred to as streaming potential, in a situation where finite temperature gradients do indeed exist. 
By way of analysing a model problem of flow through a narrow parallel plate channel, we show that 
the temperature gradients have a significant effect on the streaming potential, and, consequently, on 
the flow itself. We incorporate thermoelectric effects in our model by a full-fledged coupling among 
the electric potential, the ionic species distribution, the fluid velocity and the local fluid temperature 
fields without resorting to ad hoc simplifications. We expect this expository study to contribute 
towards more sophisticated future inquiries into practical micro-/nano-fluidic applications coupling 
thermal field focusing with electrokinetic effects. 



I. INTRODUCTION 

Streaming potential is one of the four primary elec- 
trokinetic phenomena, the other three being electroos- 
mosis, electrophoresis and sedimentation potential. The 
genesis of these phenomena is contingent on the devel- 
opment of an electrical double layer (EDL) which refers 
to the space charge distribution in a dielectric medium 
together with the electrified surface in whose immediate 
vicinity such distribution is established through a balance 
between Coulombic and entropic interactions [TJ. The 
particular phenomenon of streaming potential is, how- 
ever, set apart by the fact that its manifestation does 
not depend on the application of an external field (unlike 
electroosmosis and electrophoresis) nor does it involve 
the transport of particles bearing such electrified surfaces 
(unlike electrophoresis and sedimentation potential). As 
long as flow (actuated through simple mechanical actua- 
tion) of fluid, bearing the space charge distribution, takes 
place past the electrified surface, a streaming current will 
necessarily be manifested, and given a scope of charge 
accumulation through the specific geometry, a streaming 
potential may also be expected to be manifested together 
with such current. 

Indeed, it is because of this apparent simplicity that 
ever since the discovery of this phenomenon more than 
150 years ago by Quincke 0, streaming potential has 
been found to be a key element in the explanation of 
various phenomena in areas as diverse as physiological 
EHS] to geophysical [6j . It has also been used in a host of 
applications in the colloidal science realm; for instance, 
zeta potential measurement and electrokinetic character- 
ization of surfaces [TUllj. Most recently, it has been at 
the forefront in the world-wide research agenda of innova- 
tive energy conversion techniques [T2H32I- Nevertheless, 
despite the long and what would otherwise seem an 'es- 
tablished' history [40 - 45 , ongoing research efforts con- 
tinue to further our fundamental understanding as well 
to extend the possibilities of its use in tandem with newer 



surface and flow characteristics [461465J . 

Interestingly, such research and modeling efforts have 
been carried out, almost without exception, with the un- 
questioned assumption that the systems under consid- 
eration are isothermal. This is despite the ubiquity of 
non-isothcrmal natural settings, and also inspite of the 
fact that even in laboratory settings, perfect isothermal 
conditions are difficult to realize in practice. Such an 
approach is paradigmatic of most of the electrokinetic 
modeling efforts inasmuch as it pertains to the micro- 
and nano-fluidic context. To be sure, there do exist 
numerous works which are indeed concerned with vari- 
ations of temperature. However, these involve simplis- 
tic one-way couplings [6"6ff71"] . or, at best, couplings be- 
tween momentum and energy through temperature de- 
pendent thermophysical properties [72-74 . Importantly, 
they do not consider the fundamental dependence of the 
ionic fluxes on the temperature variation. This situ- 
ation is rather surprising particularly when considered 
in the context of the rich theory that already exists to 
model coupled transport of momentum, heat and mass 
based on general non-equilibrium thermodynamic princi- 
ples [75j [76] . Such general theories have routinely been 
adapted to represent various transport phenomena in- 
volved in membrane technology 77 . Investigations in- 
volving non-isothermal transport are also quite common 
in electrochemical systems under the purview of ther- 
moelectrochemistry |78| . particularly at high tempera- 
tures [79 . On another front, certain fundamental mod- 
elling frameworks have also been developed under the 
purview of colloidal science [801 EI] ■ Despite the obvi- 
ous commonalities that exist between these areas and 
micro and nano-fluidics, the necessity to model non- 
isothermal electrokinetic transport through these micro 
and nano-channels, in general, and streaming potential 
mediated non-isothermal transport in particular has re- 
mained largely unattended [82 . 

It is only in the last few years that this situation has 
been remedied to a great extent through a series of works 
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[8"5H9"5] where the fundamental influences of temperature 
gradient on the flux of the ionic species itself have been 
explicitly captured following the strategy by Guthrie et 
al. [M] Notable precursors to these works were the ones 
by Ruckenstein JJ5] and Morozov [95] . These works have 
also motivated the application of such gradients for con- 
trolled manipulation of particle motion in microfluidic 
channels [97]. While the modeling efforts in these works 
are directed towards the transport of colloidal particles 
and are, hence, reminiscent of the aforementioned ear- 
lier works [SHI |5T] j important differences do exist in the 
approach itself 87J. In contrast to the previous works 
which involve a dependence on the enthalpy in the veloc- 
ity expression, these recent efforts express the velocity 
in terms of the temperature gradient in an explicit and 
straightforward way. This is done by taking into account 
the contribution of the Soret effect directly in the flux. 
Pertinently, the Soret effect refers to the motion of a 
particle when placed in a temperature gradient |98l 199) . 
When the particle itself is charged (as it in the case of an 
ionic species) such transport may work in tandem with 
conventionally recognized responses to concentration and 
electrical potential gradients to generate novel flow char- 
acteristics. 

The influences of the rich interplay among these vari- 
ous factors under conditions of dynamic equilibrium are 
particularly intriguing in the case of streaming poten- 
tial mediated flow. This is because it holds the possibil- 
ity of unveiling unexpected trends in the otherwise rou- 
tine, primarily pressure-gradient actuated flow through 
the simple application of a temperature gradient. This 
motivates the primary objective of the current study: To 
investigate the influence of an externally applied temper- 
ature gradient on a pressure-gradient driven flow which 
also results in the generation of a streaming potential. 
Notwithstanding the apparent simplicity of such an ob- 
jective, the scope of this investigation is far-reaching. An 
insight into the complexities involved is in order. 

The overarching condition for setting up the dynamic 
equilibrium is that the total ionic current across the cross- 
section of the flow conduit should necessarily be zero. In 
the absence of any temperature gradient, such a condition 
results in the manifestation of a streaming potential field 
through a balance of the streaming current and the con- 
duction current. However, when a temperature gradient 
is applied, the consideration of the Soret effect may alone 
add two levels of complexity. The first rather straight- 
forward one is the motion of the particle in response to 
the temperature gradient: its direction depending on its 
thermophobicity or thermophilicity, and its intensity on 
the magnitude of the Soret diffusion coefficient which is 
effectively a function of the ionic heat of transport. The 
second level of complexity stems from the different val- 
ues of such ionic heats of transport of the cations and 
the anions which leads to a separation of charge paral- 
lel to the direction in which the temperature is applied, 
and thus sets up its own potential difference. This is re- 
ferred to as the Seebeck effect and the electric field asso- 



ciated with the potential difference is known as the ther- 
moelectric field. The most crucial thing to note is that 
under the conditions of dynamic equilibrium, a certain 
potential difference is generated along the flow direction. 
This potential difference is, however, devoid of any man- 
ifest identity of its genesis in either the Seebeck effect or 
the streaming potential phenomenon, and is solely a self- 
consistent manifestation of the combined interplay of the 
various mechanical, chemical, thermal and electrical fac- 
tors. The fluid flow that is ultimately generated is thus 
a function of the applied pressure-gradient, the applied 
temperature gradient and the self-consistently generated 
electrical potential difference. 

Notwithstanding the deep synergy among the various 
factors at work, we are able to delineate the specific ex- 
tent to which the Soret effect, in tandem with the Seebeck 
effect (that leads, in turn, to the thermoelectric field), 
influences the streaming potential mediated flow. The 
primary finding of this work is that with a sufficiently 
strong Seebeck effect, it may be possible to either aug- 
ment or negate the volumetric suppression of the primary 
pressure-driven flow due to the streaming potential field. 
This is particularly important because such reversal of 
flow nature is achieved simply on the basis of the elec- 
trolyte nature without changing the applied temperature 
gradient. 

The remaining part of this article is organized as fol- 
lows. In Section II, we describe the the model problem as 
a way of analysing the influence of temperature gradient 
on streaming potential modulated flows. We also outline 
the general equations governing the electrical potential, 
the ionic species distribution (incorporating thermoelec- 
tric effects) , the fluid-flow equations and the energy equa- 
tion that need to be solved in a coupled way for a resolu- 
tion of the intrinsic interdependence among the various 
fields. Additionally, in this section, we derive the electric 
field associated with the streaming potential incorporat- 
ing such thermoelectric effects. In Section III, we adopt 
a non-dimensional scheme and present the dimensionless 
versions of the coupled equations of Section II. In Sec- 
tion IV, we discuss the boundary conditions. In Section 
V, we report and discuss the major findings of our in- 
vestigation. Finally, in Section VI, we draw important 
conclusions based on these findings. 



II. MATHEMATICAL FORMULATION 

Pressure-driven transport of a Newtonian fluid, con- 
taining symmetric electrolytes (z + = z_ = —z), through 
a long, parallel-plate channel of height 2H, length L 
and width W (W » 2H) is considered here. The 
flow is actuated by a constant axial pressure-gradient, 
P x = —dP/dx. The fluid enters the channel at temper- 
ature To with a uniform velocity. The tip of the plate 
where the fluid enters is maintained at temperature T\. 
There is an imposed linear temperature gradient on both 
the plates. The temperature of the tip of the channel 
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FIG. 1. Schematic of the problem geometry 



where the fluid leaves is Tj. To simplify the problem 
the following assumptions are made: steady, incompress- 
ible, and laminar flow of the electrolyte is considered, 
the thermophysical properties of the fluid are considered 
to be independent of temperature variations, the ionic 
species are assumed to behave as point charges, and the 
zeta potential, £, is uniform along the channel walls. 



A. Potential distribution 

From the theory of electrostatics the potential (due to 
electrokinetic potential and thermoelectric potential), ifr, 
in the diffuse layer of EDL is governed by the Poisson 
equation which can be described as: 



Pe 

e 



(1) 



where e is the permittivity of the medium. For z : z 
symmetric electrolyte, the charge density, p e may be ex- 
pressed as: 



p e = ez (n 4 



(2) 



where e is the electronic charge magnitude, n + and n_ 
are, respectively, the number density of the positive ions 
and the negative ions. Then the Poisson equation may 
be written as: 



-) 



B. Species Transport 



(3) 



The general species transport equation routinely in- 
voked in continuum descriptions of electrokinetic phe- 
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Vj±, 



(4) 



where j± is the flux of the cation and the anion. Here, 
an important assumption is that there is no source term 
due to generation or consumption of any ionic species by 
any bulk reaction within the electrolyte. In conventional 
isothermal treatments, this flux arises as a combination of 
diffusion, electromigration and advection [lOOj. However, 
in the present case, there is an additional contribution to 
the flux from the temperature gradient. Physically, this 
comes about due to the propensity of a species to move 
in response to gradients in the temperature field, and is 
referred to as the Soret effect. This has rich implications 
in the resulting physical scenario. For, this motion of 
the ionic species due to the temperature field may itself 
result in the generation of a potential especially when 
the thermal diffusivities of the ionic species, with dif- 
ferent polarities of charge on them, are different. The 
generation of a thermoelectric potential is analogous to 
the Seebeck effect. The total electric potential which in- 
fluences the flux is, thus, a combination of the intrinsic 
electrokinetic screening potential and the thermoelectric 
potential together with the streaming potential which is 
induced parallel to the flow direction as a result of the 
streaming motion of the ions associated with the fluid 
flow . The resultant species transport is, therefore, an 
intimately coupled manifestation of this combined elec- 
tric potential, the temperature gradient, the gradients in 
the concentration (which is itself coupled to the potential 
through the Poisson equation), and the advection asso- 
ciation with the fluid motion (which, as we show later, 
is, again, influenced by the electric potential and concen- 
tration gradients). To address this highly coupled phe- 
nomenon, we first note that the flux of the ionic species 
may be expressed as [53] : 



j+ = n + u - 



j_ = ri_u — D- 



k B T 2 k B T 



k B l z KrI 



(5) 
(6) 



where u is the mean fluid velocity, k B is the Boltzmann 
constant, T is the absolute temperature, E s is the electric 
field associated with the streaming potential, D± denote 
the diffusitivities of the posistive and the negative ions, 
and Q± denote the ionic heat of transport of the posi- 
tive and the negative ions. Then, at steady state, the 
transport equations for ions may be written as: 
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o. 

= 0. 



(7) 
(8) 



C. Streaming Potential Field 

The downstream migration of ions due to the flow pri- 
marily actuated by the imposed pressure gradient gives 
rise to a current known as the streaming current (J s ). 
However, in the stationary state, this convective trans- 
port of ions sets up its own electric potential known as 
the streaming potential. The electric field (E s ) associ- 
ated with this streaming potential generates a current, 
known as the conduction current (J c ). Much like the 
electroosmotic flow situation, this conduction current, in 
turn, leads to a fluid flow opposite to the pressure-driven 
flow which was responsible for inducing the same in the 
first place. In the absence of an imposed electric field, 
the conduction current and the streaming current must 
balance each other so that the net ionic current in the 
system along the axial direction is zero; thus: 

lionic Is ~t~ -^c 0- (9) 

It is important to realize that the real physical situa- 
tion is far more involved that this rather simplistic de- 
scription might apparently indicate. This is because of 
the inricate coupling among the electric potential (with 
contributions from the electrokinetic screening potential, 



the thermoelectric potential and the streaming poten- 
tial), the concentration profiles, the temperature and the 
velocity fields. Such an intricate coupling might indeed 
include situations where locally the flow might be in a di- 
rection opposite to that expected from the primary pres- 
sure gradient drive. The only requirement that needs 
to be necessarily satisfied is that the total ionic current 
along the axial direction integrated over any cross-section 
of the channel vanishes. In order to capture this general 
requirement, we express the local ionic current in terms 
of the axial direction components of the cationic and the 
anionic fluxes: 

i = ez(j+x -j-x), (10) 

so that the total ionic current across a cross-section of 
the channel is given by: 

I-2H 

home = ez {j +x - dy. (11) 
Jo 

Thereafter, imposing the condition of the vanishing total 
ionic current, we obtain an expression for the electric 
field, E s associated with the streaming potential in the 
form: 
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(12) 



The important thing to note here is that this expression 
is not an explicit expression for E a because the velocity 
field in the first term in the numerator depends on this 
streaming potential field itself. To obtain a proper esti- 
mate of this dependence, we next move on to a descrip- 
tion of the velocity fields from the Navier-Stokes equa- 
tions representing the momentum transport of the fluid. 



D. Navier-Stokes Equation 



The axial direction velocity contributing to the flux 
and, hence, to the streaming potential field is governed by 
the x-momentum equation which under the assumption 



of low Reynolds number flow reduces to: 




-dP d 2 u d 2 u „ 

1~ + ^^2 + V^-2 + Fx - 

dx oy z ox z 



(13) 



where the total body force F x is made up of two con- 
tributions: First, Fo x due to the osmotic pressure, and, 
second, F Ex due to the Maxwell stress along the axial di- 
rection. Noting that the general expression of force due 
to osmotic pressure is: 

F = -V (nk B T) (14) 

with n = n + + n_ , the x-component of this force is: 

Fo x =-k B (n ++ n^-k B T d -^^ (15) 

It is, therefore, this contribution from the osmotic pres- 
sure which takes into account the dependence of the ve- 
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locity field on the imposed temperature gradient and the 
concentration gradients that are established as a com- 
bined consequence of the electrokinetic and thermoelec- 
tric effects. It is important to note that there does exist 
an alternative route to this explicit consideration of os- 
motic pressure contributions separate from purely hydro- 
dynamic pressure ones by considering the entire pressure 
contributions to be subsumed within the first pressure 
gradient term. Resolving the combined pressure contri- 
bution in this alternative route would, however, require 
the consideration of the transverse-direction momentum 
equation akin to semi-analytical treatments of diffusioos- 
motic phenomena (see, for example, Ref. |101L 1102] ) 
- a step which is by-passed in the present formalism. 
Next, the general expression for the contribution of the 
Maxwell's Stress to the body force is (assuming constant 



permittivity) : 



F B = eV 



(16) 



where the total potential, </>, is the sum of potentials due 
to streaming potential, <po, and the potential screening 
the surface charge, tp, so that <f> — cbo+ip. Further, noting 
that V 2 ps dy4> because 0o is constant along the trans- 
verse direction and both cf>o and tf> are assumed to weakly 
vary along the x-direction, we have = ed 2 ^ ((pa+^P)- 
Additionally, since V0o = — E s , we can write for the ax- 
ial component of this force as: 



d 2 ip 
dy 2 



dx ) ' 



(17) 



Using Eqs. 15 and |17| the x- momentum equation may 
then be written as: 



-dP 
dx 



d 2 u d 2 u 



dy 2 



dx 2 



dT 

dx 



d(n+ 



dy 2 s 



d 2 4> dijj 
dy 2 dx 



(18) 
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E. Energy Equation 

Taking into consideration the effects of axial conduc- 
tion and viscous dissipation the thermal transport equa- 
tion can be written as: 



pCpU 



dT 

dx 



= k 



d 2 T , d 2 T 



dy 2 



dx 2 



I* 



du 



(19) 



where p is the density, C p is the specific heat capacity, k 
is the thermal conductivity and u is the viscosity of the 
electrolyte. In an effort to simplify the analysis without 
sacrificing the essential physics, we consider the thermo- 
physical properties to be temperature-invariant. This 
consideration also helps us to isolate the demonstrated 
effects from electrothermal effects that may occur as a 
consequence of variation of electrical properties with tem- 
perature. 



ratio of the channel geometry is denoted by a = H/L. 
The new dimensionless variables (denoted by a tilde on 
top) are given by: 



if, 



T 



eztp 
4fc B T ' 
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(20) 



A. Dimensionless Poisson Equation 

Using the aforementioned non-dimcnsionalization 
scheme in Eq. [3j we get the dimensionless Poisson equa- 
tion as: 



d ib d ib 

or-^r- H - 

dx 2 dy 2 



= -\l< 2 {n + -n_), (21) 



III. DIMENSIONLESS FORMS 

To simplify numerical implementation, the governing 
equations are rewritten in dimensionless form using char- 
acteristic parameters of the system: the channel half- 
height H, the length of the channel L, the ionic concen- 
tration no of the bulk electrolyte, the reference temper- 
ature To, and the characteristic velocity u re f. Addition- 
ally, the potential is non-dimensionalized by the thermal 
voltage scaled by a factor of 4, ksTo/Aez identical to Ref. 
|48j . and the cationic and the anionic diffusion constants 
are assumed equal so that D+ = D- = D. The aspect 



where K = H/X and A = v / (ek B To)/(2n e 2 z 2 ) is the 
Debye length representing the characteristic thickness of 
the EDL so that K represents the penetration of the EDL 
into the bulk relative to the channel half- height. 



B. Dimensionless Species Transport Equations 

Again using the aforementioned non- 
dimensionalization scheme in Eqs. [7] and [8j the 
dimensionless transport equations for positive and 
negative ions become: 



where Pe — u re fH/D is the Peclet number, Q — 
Q + /kBT and 7 = Q-/Q+- The contribution to the ad- 
vective terms from the transverse direction velocity has 
been assumed to be negligible. 



where C — "^dpjdx re P resen ts the strength of the os- 
motic pressure relative to that of the hydrodynamic pres- 
sure. 



C. Dimensionless Streaming Potential Field 

The dimensionless form of the streaming potential field 
using the aforementioned non-dimensionalization scheme 
in Eq. |T2|is: 



-^Peh+I 2 + h + Qh 



(24) 



where the streaming potential field has been non- 
dimensionalized by ksTo/ezL, and the expressions of the 
various integrals are: 



E. Dimensionless Energy Equation 

Finally, the dimensionless energy equation by using the 
non-dimensionlization scheme in Eq. |19|is: 



aPeru 



dT 

dx 



d 2 T 
dy 2 
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, d 2 T 
dx 2 



Br R 



On 



(31) 



where Per = pC p u re fH /k is the thermal Peclet number, 
and Brji = fxu 2 e j:/kTo is a Brinkman number based on 
the reference temperature T . 
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(25) 
(26) 
(27) 
(28) 
(29) 



D. Dimensionless Navier-Stokes Equation 

The dimensionless Navier-Stokes equation using the 
aforementioned non-dimensionalization scheme in Eq. 



(181 is: 



= 2 



, a u u 
dx 2 dy 2 

8C 



C 
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K 2 



dx 

dip\ d 2 ip 
dx I dy 2 ' 



(30) 



IV. BOUNDARY CONDITIONS 
A. Wall 

The channel wall is assumed to be at a constant zeta 
potential (tp = C) with no flux of ions across it (h-j± = 0), 
with n depicting the unit vector normal to the surface. 
No-slip boundary condition is assumed: u — 0. A linear 
temperature gradient is applied along the channel wall 
(T = Ti + (T f - Tt)x/L). In dimensionless form these 
boundary conditions at the wall (y = and < x < 1) 
are expressed as: 



Q± n±dT 
k B T { f 2 dy 



= 0, 
(32) 



u=0, T = T ± {1 + x (T ratio - 1)} , 
where T ratlo = Tf/T x . 

B. Centre-line 



At the channel centre-line we assume far-stream condi- 
tion and set the electrokinetic potential to zero ("0 = 0) in 
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tandem with the electro- neutrality condition: n± = iiq. 
Exploiting the symmetry of the channel about the centre- 
line, we set the gradient of the velocity and the temper- 
ature in the transverse direction to be zero: du/dy and 
dT/dy = 0. In non-dimensional form, these boundary 
conditions at the centre-line (y — 1 and < x < 1) are 
given by: 



du 



0, n± = 1 
df 



0. 



0. 



(33) 



C. Entrance 

The fluid is assumed to enter the channel at the ref- 
erence temperature (T = To) at any arbitrary velocity 
(say, u — 0.1u re f). Assuming reservoir at the entrance 
we set ion concentration equal to the bulk concentration 
(n± = no). We also assume that the induced potential 



is zero (?/> = 0). In non dimensional form the boundary 
conditions at the entrance (0 < y < 1 and x = 0) are 
given by: 



V> = 0, 7l-| 



0.1, T=l. 



(34) 



D. Exit 

When the fluid leaves the channel it is hydrodynami- 
cally fully developed. Again, assuming the presence of a 
reservoir at the exit, we set the ion concentration to the 
bulk value, n± = uq. We also assume that the induced 
potential is zero, ip = 0. The temperature profile at the 
exit is derived from energy conservation principle (i.e. 
energy leaving the control volume is the sum of energy 
entering the control volume and the energy generated 
through viscous dissipation); this conservation principle 
is expressed in dimensional form as: 



H L 

P C p f u(0,y)T(0,y)dy-k 



dT 
dy 



H L 

dx + 

y=0 



du 
dy 



H 



dxdy = pCp I u(L,y)T(L,y)dy 



(35) 



Energy entering 



Energy generated 



Energy leaving 



In non-dimensional form, the boundary conditions at the 
exit (0 < y < 1 and x = 1): 



aPe T / u(0,y)T(0,y)dy~ 



dT 



dx + Brp> 



y=0 



du 



~ ~ du 

tp = 0, n ± = 1, ~ = 0, 
ox 



dxdy = aPe T u (l,y)T (l,y) dy. (36) 



o o 



RESULTS AND DISCUSSIONS 



The primary objective of the present work is to inves- 
tigate the influence of temperature gradients on stream- 
ing potential mediated flows. The framework required 
to model such scenarios as represented through Eqs. |21f 
36 is a weave of intrinsically coupled fields. Even with 
a number of simplifying assumptions this framework is 
manifestly more involved than those used in prior works 
which study streaming potential mediated flows in a 
number of different settings albeit, without exception, in 
isothermal conditions. Without further simplifications 
and/or approximations, this framework is not amenable 
to analytical treatments and, hence, requires the use of 



numerical techniques. For our purpose, we choose the 
finite element method as implemented in the COMSOL 
Multiphasics environment. Before setting out to report 
the results of the temperature gradient influences, we 
first establish the validity of our numerical framework by 
way of comparison of the isothermal condition results ob- 
tained from the same with those from a well-established 
semi-analytical formulation found in Ref. [JS]; perti- 
nently, this formulation has formed the basis of a num- 
ber of subsequent works dealing with streaming potential 
mediated flows [29l IMffl EH [53l EU [57l EH E3]. Next, 
we show the importance of considering the Soret effect 
by comparing results with those obtained from a consid- 
eration of temperature gradient influences solely through 
explicit manifestation of such gradients in the momentum 




FIG. 2. (Left) Numer density of counterions and coions along the cross-section at x = 0.5 corresponding to C = —0.5 and 
K — 10 and in the absence of any temperature gradient. Inset shows the corresponding electrokinetic screening potential profile. 
(Right) The corresponding axial direction velocity profile. The plots with markers represent the results from the semi-analytical 
formulation while those without represent the results from the numerical framework developed in this study. These results are 
found to be in excellent agreement. 



equation. Finally, we explicate the specifics of these Soret 
effects through their influences on the ionic species flux 
on temperature gradient mediated streaming potential 
flows. Unless otherwise mentioned, we show all results 
corresponding to constant values of £ = —0.5, K = 10, 
C = 0.828, a = 2.39 x 10~ 6 and Pe = 6.85 x 10~ 4 . 



A. Validation of simulation results 

In particular, for the sake of this validation study, we 
first switch off the influence of the temperature gradient 
in our numerical framework by setting T rat i = 1. Next, 
we restate the dimensionless axial direction velocity from 
the semi-analytical formulation of Ref. [35] in a suitably 
modified non-dimensional format to make the represen- 
tation amenable for meaningful comparison as: 



(37) 



Here the dimensionless streaming potential, E s is given 
by: 



-Pel u 



,CPe 



*2s — O-^xlhs 

where the expressions of the three integrals are: 



hs 







(n + -n_) (y 2 - 2y) dy 



hs = I {n+ + h_) dy, 
Jo 

hs = J C (n+ - n_) ( 1 - ^ I dy. 



(38) 

(39) 
(40) 
(41) 



It is to be noted that in this formulation, the electroki- 
netic potential profile is obtained by solving the Poisson 
equation: 



d 2 ^ 
dy 2 



= --K 2 (n+-n-) 



(42) 



where the ionic number densities are assumed to follow 
the Boltzmann distribution: h± — exp(=p4?/;). 

In Fig. 2 (left panel), we show the distribution of 
the number density of the counterions and the coions 
along the cross-section of the channel at x = 0.5 ob- 
tained from both the numerical implementation and the 
semi-analytical formulation. We also show, in the inset 
of Fig. 2, the comparison of the electrokinetic potential 
profile, accompanying the ion distribution, obtained from 
the two methods. Furthermore, in Fig. 2 (right panel), 
we show the velocity profiles corresponding to the distri- 
butions shown in the left panel again obtained from the 
two methods. All the results from the two methods are 
found to be in excellent agreement. This, therefore, sets 
a robust ground for using this numerical framework in 
further investigations of the influence of the temperature 
gradient. 



B. Influence of temperature gradient 

The straightforward stratagem that one might expect 
to follow while investigating any influence of the temper- 
ature would be through the incorporation of such thermal 
gradients emanating from osmotic pressure contributions 
(which otherwise remain latent in conventional isother- 
mal treatments) in the momentum equation. While such 
an expectation is not wrong it does not constitute the 
entire picture. For, even though the incorporation of 
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FIG. 3. Comparison of the numer density of counterions and 
coions considering Soret effect and that without along the 
cross-section of the channel at x = 0.5 with ( = —0.5, K — 10, 
Q — 1.388 and j = 0.153. Inset shows the corresponding 
comparison for electrokinetic screening potential profiles. 



'30.4 




FIG. 4. Comparison of the velocity profiles considering Soret 
effect and that without along the cross-section of the channel 
at x = 0.5 with (f = -0.5, K = 10. Values of Q = 1.388 
and 7 = 0.153 corresponding to a typical alkali halide |84| are 
used. These values will be subsequently referred to as Q* and 
7*- 



temperature gradients in the momentum equation and 
the consequent influences on the velocity field may in 
turn be expected to influence the ionic flux leading to 
significantly coupled manifestation of such gradients in 
the overall field distributions, it still does not take into 
account the intrinsic dependence of the ionic species flux 
on the temperature gradient. Indeed, the influence of the 
temperature gradient on the ionic flux is very deeply in- 
grained having its basis on fundamental cross-couplings 
associated with the non-equilibrium thermodynamics of 
general species transport; in particular, the Soret effect. 
This determines the 'diffusional' transport of a particle 



in response to an applied temperature gradient. 

As a first step in our discussion of the temperature 
gradient influence on the overall transport problem, we 
show, through Figs. 3 and 4, the differences in the re- 
sults obtained from the consideration of Soret effect as 
compared to the aforementioned intuitive (but physically 
incomplete) stratagem. Fig. 3 shows the distribution of 
the number density of the counterions and the coions 
across the cross-section at x = 0.5 while the inset shows 
the distribution of the potential profile across the same 
cross-section. Fig. 4 shows, again for the same cross- 
section, the velocity profile associated with the distribu- 
tions in Fig. 3. The plots with Soret effects incorporated 
are obtained for Q = 1.388 and 7 = 0.153. These values 
correspond to values of the ionic heat of transport of the 
cation, Q_|_, and that of the anion, Q_, typical for an 
alkali halide [53]; notably Q + ~ 10Q_. (In the ensuing 
discussion these values of Q and 7 will be referred to as 
Q* and 7*). While the influence of the Soret effects may 
not be evident in Fig. 3, it perceptibly clear in Fig. 4 
where it is seen to suppress the magnitude of the veloc- 
ity profile. This is particularly significant because it is 
the velocity profile which ultimately determines the vol- 
umetric flow rate and thus the throughput ratings of any 
micro- or nano-fluidic device. It is important to note that 
this suppression of the volumetric flow rate due to the 
temperature gradient influences is over and above that 
due to the streaming potential effects which invariably 
reduces the throughput in pressure-gradient driven flows 
through the generation of a back potential that drives 
a self-induced back electroosmotic flow. A discussion of 
the reason behind this suppression is in order. 

Positive values of Q and 7, originating from positive 
values of Q + and Q_, represent thermophobic nature of 
the ions. This means that within the sole purview of the 
Soret effect (independent from any other electrokinetic 
or convective influences), these ions have a tendency to 
move from the hot to the cold region. Since temperature 
increases in the direction of the pressure-gradient driven 
flow (as shown in Fig. 1), the thermophoretic movement 
(associated purely with the Soret effect) of the ions is 
in the opposite direction. In perfect anlaogy with the 
physical explanation of the phenomenon of electroosmo- 
sis where it is the electrophoretic motion of the electrical 
double layer charges which get translated into a motion 
of the fluid, so also in the current situation, the back 
thermophoretic motion of the ions results, in turn, in a 
thermoosmotic flow of the fluid. Since this temperature- 
gradient mediated back flow of the fluid is in the same 
direction as the back electroosmosis-like flow associated 
purely with the induced streaming potential, the applied 
temperature gradient is seen to suppress the fluid flow to 
a greater extent compared to the case with no considera- 
tion of the Soret effects and, consequently, no possibility 
of any such intrinsic thermophoretic motion. There is, 
however, one effect which may countervail the suppres- 
sion of the volumetric flow rate due to the Soret effect, 
and that is the thermoelectric effect. The genesis of this 
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thermoelectric effect is, essentially, the same as that of 
the streaming potential as it arises due to the accumula- 
tion of the ions which get transported in response to the 
temperature gradient. In what follows, we explicate the 
combined interplay of the Soret effect and the thermo- 
electric effect through their influences on the volumetric 
flow rate. 



0.6 
0.5 
0.4 



O 1 



0.3 



0.2 



0.1 





= 7* 




= 5 7 * 


7 


= 157* 


7 


= 207* 



0.42 
O 0.4 
0.38 



-2 

lo gio(7/7* 



-10 12 

logioCG/G*) 

FIG. 5. Variation of the dimensionless volumetric flow rate 
corresponding to the variation of Q/Q* over four orders of 
magnitude for four different values of 7/7*. Inset shows the 
variation of the dimensionless volumetric flow rate with 7/7* 
varying over four orders of magnitude. The values of £ = 
—0.5, K = 10 and T ratio = 1.17 are kept constant. 



In order to obtain a clearer picture of the combined in- 
fluences of the elect rokinetic, Soret and the as-yet unex- 
plained thermoelectric effects, we investigate the extent 
to which the values of the Soret effect parameters Q and 7 
affect the volumetric flow rate without limiting ourselves 
to the constraints of any specific electrolyte. In Fig. 5, 
we study the variation of the dimensionless volumetric 
flow rate over the half-channel cross-section (defined to 
be J Q udy) as the value of Q is varied around Q* over 
four orders of magnitude, i.e. —2 < log 10 (Q/<2*) < 2 
for different values of 7. Based solely on the explanation 
of the back flow of the ions (with the concomitant back 
flow of the fluid) presented earlier, it would be natural to 
expect that the volumetric flow rate would decrease with 
increase in the value of Q. Contrary to that expecta- 
tion, however, we observe that the value of the volumetric 
flow rate Q is practically constant (the numerical value 
actually increases slowly) until about log 10 Q/Q* ~ 0.5 
beyond which it increases fast. This increase in the volu- 
metric flow rate means that the suppressing influence of 
the Soret effect due only to the thcrmophobic back flow 
as discussed previously gets reduced. This is strongly 
indicative of the emergence of another physical factor. 
To motivate the understanding thereof, we note that in 



the case of the streaming potential the convective trans- 
port of the ions (in response to a pressure-gradient driven 
fluid flow) results in an accumulation of the ions down- 
stream leading, in turn, to the generation of a back po- 
tential which suppresses the further transport of the ions 
and the fluid, thus negating the very cause which es- 
tablishes it. In just the same way, the transport of the 
ions in response to the temperature gradient also leads 
to an accumulation of these ions albeit in the upstream 
direction. Again, just as in the case of the streaming 
potential, since the number density of the counterions is 
higher than that of the coions, there is a predominant 
accumulation of the counterions; this then leads to the 
generation of a thermoelectric field where the potential 
drops down along the direction opposite to that of the 
temperature gradient induced motion of the ions. In our 
case, this is also the direction along which the primary 
pressure gradient driven flow takes place. This thermo- 
electric field thus opposes the streaming potential field. 
Consequently, it also opposes the further flow of coun- 
terions in the upstream directin in response to both the 
Soret effect mediated thermophoresis and the streaming 
potential mediated back electroosmosis. This ultimately 
leads to a decrease in the suppression of the volumetric 
flow rate. 

A subtle point to note here is that besides the different 
number densities of the counterions and the coions, there 
is another factor which lies behind the genesis of the ther- 
moelectric field. This factor is the unequal values of the 
ionic heat of transport of the counterions and the coions 
which essentially determines the extent of their thermo- 
phobicity. A lower value of the ionic heat of transport 
of the coions than that of the counterions (in our case 
the coions are the anions and the counterions are the 
cations), represented by 7 < 1, ensures that the coun- 
terions have a greater propensity of getting accumulated 
in the upstream than the coions. It is only for the value 
7 = 1 (this is true for 7 ~ IO7*) that the thermoelectric 
field will be solely determined by the predominance in the 
number density of the counterions. The immediate de- 
duction that we can make from this discussion of the in- 
fluence of the relative values of the ionic heat of transport 
is this: as the difference in the values of the ionic heat 
of transport of the ions decreases, there will be a smaller 
manifestation of the thermoelectric effect. This is indeed 
seen to be so in Fig. 5 where the volumetric flow rate 
decreases for 7 = 57* implying that the thermoelectric 
field is now relatively weaker so that it cannot counter- 
act the Soret effect mediated volumetric flow resistance. 
Further increase in the value of 7 beyond IO7* leads to 
a reversal in the nature of the thermoelectric field. In- 
deed, with 7 now being greater 1, it is the coions which 
have a greater propensity to flow upstream compared to 
the counterions: this leads to an inversion in the polarity 
of the thermoelectric field. Consequently, this thermo- 
electric field now aids the streaming potential mediated 
back electroosmotic flow leading to a significant decrease 
in the overall volumetric flow as seen in Fig. 5 for the 
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plots corresponding to 7 = 157* and 7 = 2O7*. More- 
over, stronger the thcrmophobicity of the ions indicated 
by higher (positive) values of the ionic heats of transport, 
stronger is the Soret effect induced back flow and higher 
is the magnitude of the thermoelectric field: this explains 
the strong discrepancies in the values of the overall vol- 
umetric flow rate for higher values of log 10 Q/Q*. The 
inset of Fig. 5 clearly shows the decreasing trend in the 
volumetric flow rate with increasing value of 7 relative to 
7* for a constant value of Q — Q* . An extremely point 
to remember is that no matter what the values of Q and 
7 are, the combined consequences of the thermoelectric 
field, the Soret effect and the streaming potential medi- 
ated flows are fundamentally determined by the criterion 
that the net ionic current across any cross-section of the 
channel should necessarily be zero in the stationary state. 

VI. SUMMARY AND CONCLUSIONS 

In this study, we have presented the as- yet unaddressed 
consequences of temperature gradient on a streaming po- 
tential mediated pressure-gradient driven flow. We have 
incorporated the explicit dependence of this temperature 
gradient on the flux of the ions by considering the Soret 
effect in our modeling framework. To highlight the con- 
sequences of such consideration, we have first shown the 
differences between the cases with and without Soret ef- 
fects. It is revealed that inclusion of Soret effects does 
indeed influence the velocity profiles quantitatively. Ad- 
ditionally, we have also studied the influence of the See- 
beck effect which arises due to the differences in the ionic 
heats of transport of the cations and the anions. As the 
most important finding of our study, it is revealed that 
depending on the strength of the Seebcck effect and the 



concomitant thermoelectric field it generates due to sep- 
aration of charges, the suppression of the volumetric flow 
rate due to the streaming potential field may be aided 
or opposed. Importantly, such an effect is manifested for 
an unchanged externally applied temperature gradient, 
and is, thus, a sole manifestation of the nature of the 
electrolyte. 

These results have far-reaching practical consequences. 
Based on our findings, one may devise strategies in 
which temperature gradients may be employed for tun- 
ing the flow by aiding or opposing the streaming poten- 
tial, depending on the relative dominance of Soret effect 
and thermoelectric effect. Since our studies are based 
on temperature-independent thermo-physical properties, 
the distinction between the new effects brought out in 
this work and electrothermal effects routinely addressed 
in microfluidics literature are clear. Moving further for- 
ward, one may exploit such flow control mechanisms 
for designing separation platforms based on a compet- 
ing environment of thermodiffusion and electrophoresis. 
As delineated, these individual effects are not linearly 
super-imposable, as attributable to an intricate coupling 
between the thermal, potential, solutal, and velocity 
fields. The concerned implications, being significantly 
non-intuitive in nature, may hold the promises of ad- 
dressing a new paradigm of microfluidic devices that re- 
lies on a strong thermo-electrical coupling with the aid 
of an imposed temperature gradient, without necessarily 
appealing to electrothermal principles. Our work further 
demonstrates that this subtle conjecture may be verified 
in complete absence of any externally applied electrical 
field, by appealing to the establishment of streaming po- 
tential fields in the presence of a simple pressure-driven 
flow. 
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